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We reveal that Hopf solitons can be stabilized in rotating atomic Bose-Einstein condensates. The 
Hopfion is a matter-wave vortex complex which carries two independent winding numbers. Such 
a topological solitonic structure results from a superfluid flow of atoms simultaneously quantized 
in poloidal and toroidal directions. In the framework of a dissipative mean-field model we observe 
different unstable evolution scenarios of the Hopfions. We demonstrate energetic and dynamical 
stability of the Hopf solitons under experimentally feasible conditions. 
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I. INTRODUCTION 

Quantization of vorticity is a remarkable manifesta¬ 
tion of superfluid properties of Bose-Einstein condensates 
(BECs). Different kinds of vortex structures have been 
theoretically predicted and observed experimentally in 
atomic BECs (see e.g. and references therein): 

single vortex lines, vortex-antivortex pairs, vortex ar¬ 
rays, vortex rings, vortex knots, dimerons [HI, three- 
dimensional vortex solitons [HI, and solitary waves mov¬ 
ing along a straight vortex line Q (which are similar to 
“hoop” structures known in field theory [H). 

It is of special interest to construct stable vortex- 
soliton states with complex structure, such as vortex 
knots liliii, three-dimensional Skyrmions, [iMl and 
Hopfions. A Hopfion (or Hopf soliton) is a topologi¬ 
cal soliton with two independent winding numbers: the 
first, S, characterizes a horizontal circular vortex em¬ 
bedded into a three-dimensional soliton; and the sec¬ 
ond, m, corresponds to vorticity around the axis, per¬ 
pendicular to this circle. Hopf solitons appear in many 
fields, including field theory, optics, ferromagnets, and 
semi- and superconductors. Hopfions have been pre¬ 
dicted theoretically both in systems with multicompo¬ 
nent wave functions [THI and in single-component BEC 
0- Multicharged (m > 1, S > 1) vortex structures 
have been demonstrated [I^ to be unstable in opti¬ 
cal media. Stationary Hopfion solutions with m = 1 , 
S' = 1 of the single-component Gross-Pitaevskii-Equation 
(GPE) have been obtained numerically in Ref. [ly for 
a spherically-symmetric trapping potential. To the best 
of our knowledge, neither a dynamical nor an energetic 
stability analysis has been performed previously for Hop¬ 
fions in trapped BEGs. Very recently 0 Hopfions have 
been demonstrated to be stable in a single-component 
BEG without linear trapping potential but with a re¬ 
pulsive non-linearity, which grows fast enough from the 
center to the periphery of the trap. 

In our recent work [^ we suggested an experimen¬ 
tally feasible trapping configuration that can be used to 
create, stabilize, and manipulate a vortex ring in a con¬ 
trollable and nondestructive manner. In this paper we 
use a similar trapping potential and a rotating conden¬ 


sate which allows the formation of a stable Hopfion. We 
demonstrate both energetic and dynamical stability of 
the Hopf soliton for realistic experimental parameters. 


II. MODEL 


The properties of an atomic BEG close to thermody¬ 
namic equilibrium can be accurately described in mean- 
field approximation by a dissipative Gross-Pitaevskii 
equation (DGPE), which takes the following form in a 
rotating frame [2ll-l^ , 


{* - 7 )^ 


5 ^'(r, t) 
dt 


-^V2 + Kxt(r)-^^L, 

+g\'^{r,t)\'^ - (1) 


Here, 7 ^ 1 is a phenomenological dissipation parameter, 
O the rotation rate, is the Laplace operator, Lz = 
—ih{xdy — ydx) is the angular momentum in 2:-direction, 
g = Airh/as/M is the coupling strength, y is the chemical 
potential of the condensate, M is the mass of the atom, 
and as is the s-wave scattering length. We assume that 
the wave function is normalized to the number of atoms 



In the absence of dissipation (7 = 0), the particle number 
N as well as the total energy 

E = J ||^|V^|2 + 14xt(r)|^|2 

+ ( 3 ) 

are integrals of motion. 

In our model the external trapping potential 


V,^t{r)=Vh{r,z) + Vk{r,z) ( 4 ) 


is created by an oblate harmonic trap 

Vh{r,z) = ^M{ojy+ojy) 
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and two blue-detuned laser beams: a radial Laguerre- 
Gaussian beam LGq and an elliptic highly anisotropic 
“sheet” beam creating a tight repulsive potential in 
dimension: 

Vb(r, z) = Vq , 


where r = is the effective width of the 

sheet-beam, is the radial coordinate of the trap min¬ 
imum, and q is the topological charge of the red-detuned 
LGq optical beam. 

Let us consider a BEG cloud with N = 10^ of ^^Rb 
atoms {as = 5.77 nm) in an oblate harmonic trap with 
ujr = 27r X 100 Hz, co’;^ = 2 7r x 300 Hz , and a characteristic 
oscillatory length = ^/h/{Muj^ = 2.7 jam. The tube 
beam has parameters q = Rt = 13.5 /im, and Vq/h = 
3.4 kHz, the sheet beam has parameters Zt = 3.38/rm, 
and Vz/h = 1.8 kHz, where h is the Planck’s constant. 
This trap geometry is similar to that used in Ref. [2^ 
modified slightly to increase the stability of the vortex 
line by squeezing the trap in the 2: direction and shifting 
the toroidal void position farther from the center. Such a 
trap configuration is supposed to naturally accommodate 
the complex core structure of a Hopfion (see Fig. [1]). 



FIG. 1: (Color online) Schematic illustration of a vortex Hop¬ 
fion in a trapped BEC. We show a condensate isodensity sur¬ 
face consisting of the outer surface (green) and the interior 
toroidal void (light blue). A Hopfion is formed by two types of 
phase singularities: the vertical vortex line (blue) surrounded 
by the vortex ring (red). The white line with arrows shows 
schematically the particle flow. 


In what follows we use dimensionless units for the spa¬ 
tial coordinates (x, y^z) (x//^, ^//^, ^//^), for time t 
ujrt^ and for other quantities T ^ ^ E/{hwr)^ 

^ ^ C/lrj and Q Q/uJr. The DGPE then becomes 






(5) 


where g = Anasllr = 2.76x10 The trapping potential 
in dimensionless units takes the form 




v.e 


with the parameters a = Iz/Zt = 0.8, P = IzIRt = 0-2, 
Vq = 34, Vz = 18, A = UJzloJr = 3. 


For dynamical calculations the phenomenological dis¬ 
sipation parameter 7 = 0.02 was chosen to be constant, 
though general results do not qualitatively depend on this 
value. 


III. ENERGETIC STABILITY OF VORTEX 
HOPFIONS 

If we aim to study the dependence of the total energy 
(|3]) on the rotation rate Q then it is convenient to split it 
into two parts 


E{n) = Eo- nc. ( 6 ) 


The first part Eq is the energy in the absence of rotation 



dr. 


and the second part jC depends on the average of the 
angular momentum 

c = J 

We obtain the BEG ground-state Tg.s. by imaginary 
time propagation [2^ from Eq. (j5j) with 7 = 0 and then 
compute its energy Eqs from Eq. Then, in order 
to model a Hopfion state characterized by topological 
charges S and m, we imprint on the ground state an 
S'-charged vortex ring with coordinates (vr^Zr) and an 
m-charged vortex line shifted by ri off the center. 


^(r) = {/.(r)e*^7^ {/,(r)e*^'j"” ^g.^r) (7) 

with tand^ = {z — Zr)l{r — r^), tand/ = yl{x — n), r = 
/r(r) = tanh 

and /;(r) = tanh {\/(x — riY + 2/^/C)- The width of the 
vortex core is estimated by the healing length ^ at the 
peak density of the condensate. 

If the vortex line is straight and parallel to the z-axis 
and the vortex ring is parallel to the (x,^)-plane, the 
superfiows of the vortex ring and line are moving in per¬ 
pendicular directions. Additionally, the density pertur¬ 
bations caused by the vortex cores are tightly localized 
and do not overlap in the regions of interest. Therefore, 
the contributions to the total energy from the vortex ring 
and vortex line are independent and can be analysed sep¬ 
arately. 


E{^) = Eo^r{rr,Zr) + Eo,i{ri) - 0£i(n), 

and we take into account that the vortex ring in the 
(x,^)—plane does not contribute to the angular momen¬ 
tum in ^-direction. 

The energy analysis of the vortex ring is performed 
by putting S = 1 and m = 0 in Q- Then the energy 
E = Eo,r(^r5^r) depends only on and Zr (see Fig. |2J). 
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The vortex ring nucleation energy shows a minimum lo¬ 
calized close to the density minimum at the crossing of 
the “tube” and “sheet” beams. It is worth noticing that 
the existence of an energy minimum does not necessarily 
correlate with the existence of a density minimum in the 
ground state. Figure [3] shows that for certain intensities 
of the “tube” beam an off-center minimum exists in the 
vortex ring nucleation energy but not in the ground-state 
density. 



FIG. 2: (Color online) Density distribution of a trapped BEG 
cloud with contour lines of the vortex ring nucleation energy. 
Locations of energy extrema are denoted with “plus” (max¬ 
ima) and “minus” (minimum) signs. 




FIG. 3: (Color online) Absolute value of the wavefunction 
(left) and vortex ring nucleation energy (right) for different 
values of Vq {vq = 25 — blue solid line, Vq = 34 — red dashed 
line, Vq = 40 — green dotted line). 


The energy associated with the vortex line is obtained 
by putting S = 0 and m = 1 in Eq. ©• This energy 
E{Q) = Eo^i{ri) — ClCi{ri) depends only on Q and ri. In 
analogy to the analysis of Fetter we can define two 
characteristic rotation speeds. The first, is chosen 
such that E = F^g.g. at n = 0. It means that for Q > 
Qc a state with the vortex at the origin is energetically 
preferable over a state without vortex. This rotation rate 


can be defined as 




Eo^i{0) - F^g.s. 

A(o) ■ 


Note that Ci{0)/N = m = 1. Therefore, in order to 
obtain this critical velocity only £’o,^(0) has to be calcu¬ 
lated. 

The second characteristic rotation rate ftrn corre¬ 
sponds to a vanishing curvature of E{ri) at r/ =0: 


— 


d^Eo^iin) /d^Cin) 


drf 


drf 


ri^O 


This rotation rate defines the region Q > Clm where the 
energy shows a minimum at r/ = 0 and the vortex line 
at the origin can be stable (or metastable). Due to our 
specific trap configuration there might be another energy 
minimum for the vortex line which is close to the peak 
intensity of the “tube” beam. But this minimum is not 
relevant for the stability of a Hopfion as it correlates with 
a stable position of the vortex ring. As the line and the 
ring have to be spatially separated, we demand the exis¬ 
tence of a minimum at the origin as a situation favourable 
for a stable Hopfion. 

As is seen from Fig. 01 in the suggested trap configura¬ 
tion, a vortex state can be energetically preferable with¬ 
out being stable at the center of the trap. This means 
that if we need to stabilize the vortex core at the axis, 
the rotation rate has to be considerably higher than Qc- 
Therefore Clrn > which is opposite to the case of a 
spherical condensate studied in Ref. . As a result we 
consider the obtained value = 0.28 as a minimal ro¬ 
tation rate at which the vortex line can be stable at the 
center of the trap and, therefore, at which we can also 
expect the overall stability of the Hopfion complex. 



FIG. 4: (Golor online) Nucleation energy of the vortex line 
depending on the displacement from the trap center for dif¬ 
ferent values of the rotation rate. 
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IV. DYNAMICS OF HOPFIONS 

For the solution of the time-dependent GPE we use 
a split-step Fourier transform method jUj. The initial 
state for the simulations takes the same form Eq. © as 
in the energetic stability analysis of the previous section, 
i.e. a stationary state with imprinted vortex ring and line 
(S' = 1, m = 1). Their initial positions are considerably 
displaced from their equilibrium (by up to O.blr). We 
study the dynamical stability of this Hopfion solution for 
a wide range of rotation rates ft. 

We observe two different instability scenarios of the 
Hopfion dynamics. If the rotation rate is low or the ini¬ 
tial position of the vortex line is too far away from equi¬ 
librium, the system is unstable with respect to reconnec¬ 
tion of the line and ring cores. As a consequence, several 
line vortices are created which, in turn, can be stable 
and form a lattice at the off-center energy minimum (see 

Fig. El). 




FIG. 5: Snapshots of the unstable Hopfion dynamics for the 
rotation rate Q = 0.2. In the initial state S' = 1, m = 1. 
Reconnection happens at t = 0.135 s. The four vortex lines 
in the final state are stable. 

The second instability scenario is observed at high ro¬ 
tation rates and is independent of the initial displace¬ 
ment of the vortices. Due to a low energy barrier (see 
Fig.lH), additional vortices can be nucleated at the outer 
edge of the BEG cloud and then penetrate into the con¬ 
densate destroying the ring vortex and again forming a 
regular vortex chain at the off-center energy minimum 
(see Fig. Ej)- 

It is however possible to choose a rotation rate that 
is high enough for stabilization of the vortex line at the 
center but at the same time low enough that additional 
vortices are not generated at the outer edge of the con¬ 
densate. Then the vortex Hopfion may be dynamically 
stable for experimentally traceable time intervals which 
are comparable to the life-time of the condensate (see 
Fig. [7]). For a strong initial displacement of the vor¬ 
tex line of 0.5 /r we found a region of stability between 
0.225 < ft < 0.245. It is worth noticing that the vor¬ 
tex line relaxes to a position slightly off the center of the 




FIG. 6: Snapshots of the unstable Hopfion dynamics for the 
rotation rate D = 0.3. In the initial state S' = 1, m = 1. 
The Hopfion complex itself is stable, but starting from t ~ 
0.46 s vortices from the outside start to penetrate the cloud 
eventually destroying the vortex ring. 

trap even though the energy minimum is located there. 
However, the energetic stability analysis does not take 
into account bending of the vortex cores, which are ener¬ 
getically preferred over straight vortices. This may also 
explain the fact that the vortex line is dynamically sta¬ 
bilized at a rotation rate lower than predicted by the 
energetic analysis. 




FIG. 7: Snapshots of the stable Hopfion dynamics for the 
rotation rate D = 0.24. In the initial state S — m — 1. 
Shapes of the vortex cores persist through all the time of 
simulation. 

It is also instructive to analyze the dynamics of Hop- 
fions with a topological charge of the vortex line m > 1. 
As multiply charged vortices in BEG are essentially un¬ 
stable, the vortex line splits into a tangled complex of 
several singly charged vortices (see Fig. 0). Because these 
vortices repel each other they can not all be simultane¬ 
ously stabilized at the center of the trap. They rather 
relax to the energy minimum at the peak intensity of the 
“tube” beam. This inevitably leads to reconnection with 
the vortex ring and the observed dynamical scenario is 
similar to the instability scenario of singly charged vortex 
lines at low rotation rates (compare Figs. El and [8|). 

Similarly, the Hopfion with a multicharged vortex ring 
{S > 1) can not be stabilized in a suggested trap geome¬ 
try (see Fig. 0. The multi-charged vortex ring splits into 
two leapfrogging singly-charged rings, that spiral around 
the energy minimum drifting away from it until they are 
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FIG. 8: Snapshots of the dynamics of the Hopfion with 
doubly-charged vortex line (S' = 1, m = 2) for the rotation 
rate Q = 0.25. 


destroyed by growing Kelvin waves. It is remarkable that 
in sharp contrast to studies of multicharged (S > 1) Hop- 
fions in optical media [TsI we do not observe the forma¬ 
tion of vortex knots in our simulations of the dissipative 
dynamics of Hopf solitons in BEG. 




FIG. 9: Snapshots of the dynamics of the Hopfion with 
doubly-charged vortex ring (S = 2, ?n=l) for the rotation 
rate Q = 0.25. 


The destruction of the Hopfion in the proposed setup 
always leads to formation of a vortex chain along the off- 
center energy minimum. The number of vortices in this 
chain is defined not only by the rotation rate but also by 
the total energy released by destruction of the Hopfion 
complex (note a six-vortex final state in Fig. [9]). 

The main aim of the present work is a demonstration 
that Hopfions can be stabilized in experimentally feasi¬ 
ble one-component trapped BECs. Suitable parameters 
sets for the stabilization of Hopfions depend on the spe¬ 
cific experimental realization of BEG. We observed that 
Hopfions are more stable in highly oblate traps since vor¬ 


tex line bending and tilting are suppressed under these 
conditions [26|. Eurthermore, imprinting the vortex line 
closer to the axis extends the stability region towards 
lower rotation rates. 


V. CONCLUSIONS 

We show that Hopf solitons can be energetically and 
dynamically stable in a rotating trapped atomic BEC. In 
the framework of a dissipative mean-field model we inves¬ 
tigate different dynamical regimes of Hopfion stability. It 
turns out that within a certain interval of BEC rotation 
velocities a Hopfion may be stable. Eor rotation rates 
lower than this interval the vortex line and ring deform 
and finally reconnect. Eor over-critical angular velocities, 
vortex lines are nucleated at the outer edge of the conden¬ 
sate and then reconnect with the vortex ring. Multiply 
charged vortex lines split into singly charged vortex lines, 
which then reconnect with the ring component. 

Decay of the Hopfion leads to the formation of a vor¬ 
tex lattice. The kinetic energy of the ring component is 
transferred via reconnection to the vortex lattice. As is 
well known, the experimental detection of vortex rings is 
much more challenging than observation of the holes in 
the atomic cloud produced by vortex lines. By counting 
the number of vortex lines in the vortex lattice one ob¬ 
tains clear experimental evidence of the presence of the 
vortex ring component. 

We hope that the results presented in this paper will 
stimulate experiments to observe Hopf solitons in atomic 
BECs. Moreover, investigation of the Hopfion in the well 
controlled environment of atomic physics could help to 
elucidate the properties of these fascinating topological 
structures in other fields including superconductors, fer- 
romagnets, and quantum field theory. 
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